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Developed Navier-Stokes turbulence is simulated with 
varying wavevector mode reductions. The flatness and the 
skewness of the velocity derivative depend on the degree of 
mode reduction. They show a crossover towards the value of 
the full numerical simulation when the viscous subrange starts 
to be resolved. The intermittency corrections of the scaling 
exponents £ p of the p th order velocity structure functions seem 
to depend mainly on the proper resolution of the inertial sub- 
range. Universal scaling properties (i.e., independent of the 
degree of mode reduction) are found for the relative scaling 
exponents p p , q = « p - (3p/3)/(( q - C3<?/3). 

Even today fully developed turbulence is hard to ac- 
cess through full numerical simulations of the Navier- 
Stokes equations because the number of degrees of free- 
dom increases with the Taylor-Reynolds number roughly 

as Re 9 / 2 @Ql • Consequently, models and approximations 
of the Navier-Stokes dynamics with a reduced number of 
degrees of freedom are considered. Models which em- 
body the cascade type structure of turbulence enjoyed 
increasing popularity in recent years, e.g., the so called 
GOY model Closer to the Navier-Stokes dynam- 

ics is its reduced wavevector set approximation (REWA) 
g-§. REWA uses a reduced, geometrically scaling sub- 
set of wavevectors on which the Navier-Stokes equation 
is solved. Very high Taylor-Reynolds numbers up to 
Re x = 7 ■ 10 4 Rl can be achieved. 



However, a priori it is not clear whether these models 
and approximations are in the same universality class as 
the Navier-Stokes dynamics itself, as small scale struc- 
tures corresponding to the high k modes are not be fully 
resolved. If inertial subrange (ISR) scaling properties de- 
pend on details of the viscous subrange (VSR) as specu- 
lated for the GOY model , a cascade type approach 
towards fully developed turbulence may not give the cor- 
rect inertial range scaling properties. Moreover - and as 
we will see more importantly - in these models the phase 
space has a different representation than in 3D Navier- 
Stokes turbulence. Indeed, detailed REWA calculations 
j^U for the scaling exponents ( p of the p th order longi- 
tudinal velocity structure functions 

•> (p) (r) = ((ui(x + ra) - Ui(x)) p ) cx r Cp , 
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i = 1,2,3, show much smaller (but non vanishing [[7|,|E)|) 
deviations SC P = Cp ~ P/3 from their classical values 
Cp = p/3 ( "K41" ) than those from experimental measure- 
ments [ ^2p^ ] or full numerical simulations (for Reynolds 
numbers up to Re\ w 210) [ld|JT5[|. Also, the flatness 
Fi = ((ftu,) 4 ) / '((diUi) 2 ) 2 is w 3.15§ for all Re x in con- 
trast to experiments and full simulations where it seems 



to increase |T(| with Re\ . Analogous results hold for the 
skewness Si = ((diUi) 3 ) / ((diUi) 2 ) 3 ^ 2 . On the other hand, 
REWA may well represent the "correct" large Re\ > 10 3 
limit where F and S are speculated to become indepen- 
dent of Re x 0]. 

In this letter we systematically analyse how the scaling 
properties change with an increasing degree of wavevec- 
tor mode reduction, i.e., we examine the transition from 
full numerical simulations to reduced wavevector set ap- 
proximations. Since full simulations are possible only for 
low Re\ values, the present calculations are restricted 
accordingly, even though REWA was constructed for the 
large Re\ limit. There is at most a short ISR. However, 
the extended self-similarity method (ESS) jl3| allows us 
to extract scaling exponents. 

The aim of the work is to better understand the ori- 
gin of intermittency scaling corrections. Two views are 
discussed: The meanwhile classical multifractal picture 
(see e.g. ref. for a review) in which intermittent fluc- 
tuations build up in the ISR and the Leveque-She re- 
flection picture |T1J in which ISR quantities depend on 
VSR properties. The importance of the latter mecha- 
nism could be shown for the GOY model |ll|]lOf| . How- 
ever, for 3D Navier-Stokes turbulence with different kind 
of hyperviscosity jl5) , no dependence of the ISR scaling 
exponents £ p on the kind of hyperviscosity could be de- 
tected. Our analysis seems to support this result. VSR 
effects on £ p could not be identified. Our interpretation 
is that the proper local phase space resolution is of prime 
importance for the correct representation of the scaling 
corrections 5£ p . 

We now describe our analysis in detail. The 3D incom- 
pressible Navier-Stokes equations are numerically solved 
on a N 3 grid with periodic boundary conditions. Spher- 
ical truncation is used to reduce aliasing. We force the 
system on the largest scale (wavevectors k = (0,0, 1)/L 
and permutations thereof) as e.g. described in ref. M. 
Units are fixed by picking the length scale L = 1 and 
the average energy input rate (= the energy dissipation 
rate) e = 1. The Taylor- Reynolds number is defined as 
Re x = Ui, rms X/p, where A = u\ irms /{d\Ux) rms is the 
Taylor length and v the viscosity. Our results refer to 
N = 60 and v — 0.009, corresponding to a resolution 
of scales r > 2nL/N w 3.677 and Re x ~ 100. Time 
integrations of about 60 large eddy turnover times are 
performed. Averages are taken over space and time. We 
also did shorter runs for N = 80 and longer runs for 
N = 48 which gave the same results. 

As our key parameter we now introduce the wavenum- 
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ber ks with 1 < ks < k rnax = N/2, characterizing the 
degree of mode reduction: For a simulation with given 
ks all wavevectors with |fc| < ks and scaled replica 2 l k, 
I = 1,2,3,..., thereof are considered; the mode ampli- 
tudes of the remaining wavevectors are put to zero. The 
choice ks = k max — N/2 corresponds to a full simula- 



tion, ks ~ 2 is our former REWA calculation |5|-|8[. For 
those calculations a pure spectral code could be used; 
here, because of the huge increase of couplings, a pseu- 
dospectral code as described in [p~8|Jl4| was employed. 




FIG. 1. Skewness S and flatness F as functions of ks- 
The error bars express the statistical differences of the values 
Si and Fi for the three space directions i = 1,2, 3. With the 
chosen lattice resolution 5* and F roughly reach their satura- 
tion value. The crosses on the very left refer to the REWA 
calculation with fcs = 3, but with a full VSR resolution for 
k > 9. 

Figure [I] shows the skewness S and the flatness F as 
a function of ks- A crossover at ko ~ 1/(477) ~ 9 can- 
be identified. Here, ku denotes the wave number with 
maximal dissipation rate, where massive viscous damp- 
ing starts in the spectrum. For fc^ < kp the flatness 
and the skewness essentially remain on their REWA val- 
ues. But at ks > ko they start to drastically increase 
towards their saturated values corresponding to the full 
simulation. 

Figure ^| shows the compensated structure function 
D^\r)/[D*^\r)] 2 vs D*^\r). This kind of plot allows 
for a better detection of local deviations from scaling than 
the standard ESS @ plot vs D*^ . We find that for 
ks > 5 the value <5Ce ~ —0.22 is always a good fit in the 
large r regime between 2tt /ks and L. This scaling regime 
shrinks for decreasing ks and vanishes below ks ~ 5 as 
then 2-K/ks essentially collapses with the external length 
scale L. 

Figure |^ suggests that at least for small Re\ for 
10?7 <C 2-n/kB <C L (a condition which never is re- 
ally reached in our small Rex simulations; the simula- 
tion for /cb = 6 is closest to it, see in particular figure 
^b) there are three ranges: The (underresolved) VSR 
r < IO77 where of course cx (L>( 3 )) 2 , a REWA ISR 
in the underresolved regime [IO77, 2tt /ks] with very small 



but nonvanishing (note the nonzero slope in figure |2Jb in 
that regime) intermittency corrections [Q,^) , and the fully 
mode resolved Navier-Stokes ISR [2ir/kB, L] with the in- 
termittency corrections S(^q = —0.22 as in full numerical 
simulations. This prompts the conclusion that it is the 
local phase space resolution and not a proper VSR reso- 
lution which is essential for the correct representation of 
scaling corrections. 
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FIG. 2. (a) Com- 

pensated ESS type plot D (6) (r)/(D* (3) (r)) 2 vs D* (3) (r) for 
various /cb. The slopes in this type of figure are the <5£ p 's; a 
horizontal line thus means K41. The open arrows point to- 
wards r — IO77 and the filled ones to r = 27r/fes. The outer 
length scale L is the very last point of each curve. Figure (b) 
shows an enlargement of (a) which illustratively shows the 
scaling with £g = —0.22 in the [27r/fc{,, L] regime. 

To further support this statement, we performed a sim- 
ulation with ks = 3, but in addition a full resolution of 
all modes k > 9, i.e., of the VSR. This curve is also shown 
in figure ^. Indeed, there are hardly any scaling correc- 
tions in the ISR, <5Ce ~ 0. On the other hand, as expected 
both the flatness and the skewness are much bigger than 
in the REWA calculations || as now the VSR is bet- 
ter resolved. We added the corresponding data points in 
figure 0. 

As is well known, the Navier-Stokes intermittency cor- 
rections are well fitted by the She-Leveque (SL) model 
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^o(|(l-/3 3 )-(l-/3 p )) 



(2) 



with the parameters (3 and Co, which in ref. |19[ where 
suggested to be (3 = (2/3) V 3 and C = 2. In eq. §) 
we already used the restriction £3 = 1 to eliminate a 
third parameter which was introduced in the original 
work Jll],^l| • The parameter Cq was related to the rate at 
which the probability to find the most intermittent events 
decays in the large k limit and also interpreted as the the 
codimension of the dissipative structures |l9|,^l],[llj . If in 
3D Navier-Stokes turbulence these are ID filaments, we 
have Co — 2. This interpretation also works for REWA 
turbulence [ks = 3) The dissipative structures are 

nearly 3D because of the lack of large wavevector resolu- 
tion. Therefore, Co ~ and according to (0) £ p « p/3, 
in agreement with the numerical results [Qj9|. However, 
the interpretation of Co as codimension of the dissipative 
structures seems to be at variance with the simulation 
with IO77 <C 27r/fce <C L, where we find the 3D Navier 
Stokes values for 5£ p in the [2ir/kB, L] regime in spite of 
the poor VSR resolution which determines the dimension 
of the dissipative structures. 

Whereas the scaling corrections 8£ p do depend on the 
local phase space resolution, their ratios 



Pp,q 



C P - C3P/3 

( q - C39/3 



(3) 



do not. These exponents are the relative scaling exponent 
of two compensated structure functions f22fj 



G {p) {r) = D ip \r)/(D* {3 \r)) p / 3 . 



(4) 



Indeed, the scaling of vs is perfect in the whole 
range we resolve, as seen from the so called general- 
ized ESS (=GESS, ||) figure |. The reason is that in 
the GESS plot all data points of the VSR collapse since 
q(p) ( r ) i s constant for r < IO77. The p p . q exponents could 
be obtained from plots like that in figure ||. Alternatively, 
they can be obtained via the 5( p from straight line fits to 
compensated ESS type plots as in figure which turns 
out to result in a smaller statistical error. Though ob- 
viously the 5( p in figure ^ depend on the chosen D*( 3 ) 
range, their ratios p Piq do not depend on this range. 

From figure || we see that the also do not depend 
on ks- Similar universality has recently been found by 
Benzi et al. [^3| for the GOY shell model where the ( p 
depend on the spacing between the shells, but the p Ptq 
do not. 

The relative scaling exponents p p%q are rather well de- 
scribed both by the SL model, predicting 



Pp,q — 



(l-^)-( P /3)(l-/3 3 ) 
(l-/J«)-(g/3)(l-/33) ! 



(5) 



which only depends on (3 and not on the codimension of 
the dissipative structures Co any more, and by the log 
normal model (LN, see ref. M for a detailed discription), 
predicting 




FIG. 3. GESS type plot G (6) (r) vs G (4) (r) for various 
ks- The inset focuses on the REWA scaling range for much 
larger Re\ rss 7 ■ 10 4 Q. The G-functions in the inset are 
calculated from the total structure functions, not from the 
longitudinal ones as in the main figure. Even for this huge 
Reynolds number the scaling range is very short because of 
the small intermittency corrections 5( p . Note that in this type 
of plot the whole VSR collapses into the upper right point of 
the curves. 
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Our values are roughly 6% below those from the LN 
model and 6% above the SL prediction with the origi- 
nally suggested — (2/3) 1 / 3 . By taking the suggestion 
of Chen and Cao pfj| j, j3 — (7/9) 1 / 3 , the agreement can 
be improved. 

Our results for both the skewness S, the flatness F, 
and for the scaling corrections 5£ p are consistent with 
Cao et al.'s jl5| full numerical Navier-Stokes simulations 
for normal and hyperviscous damping (— ) h vV 2h u (h = l 
means normal viscosity). Our large wavevector reduction 
is similar to a kind of hyperviscosity. Indeed, the larger 
h the smaller the flatness is. As in the calculations pre- 
sented here, the ISR scaling properties S( p are hardly 
affected by the kind of hypervisocity. 

In summary, we repeat that in 3D Navier-Stokes turbu- 
lence the main origin of intermittency corrections seems 
to be the proper resolution of the phase space at the 
scale of interest. Reflections from the VSR seem to be 
of minor importance. How come that the energy flux re- 
flected from the VSR |ll|,^0[ reaches so far in the ISR for 
the GOY model but apparently not for 3D Navier-Stokes 
turbulence? We speculate that in 3D the phases of the 
modes are subjected to far more fluctuations than in the 
ID GOY model. Therefore, coherences get destroyed eas- 
ier. Some coherence however must remain, otherwise no 
energy could be transported downscale. 
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FIG. 4. Relative scaling exponent p§ t 4 as a function of 
our mode reduction parameter fcs. Here, r\ is the Kolmogorov 
inner scale. The error bars take care of both the quality of 
the linear regression of the scaling laws and the deviations 
from isotropy. The triangle on the very left refers to the 
REWA calculation from ref. g for much larger Re\. The 
dotted line shows the SL prediction |l9| from equation (^) 
with /3 = (2/3) 1/3 » 0.8736, i.e., p e>4 = 4.14, the dashed one 
the prediction from the LN model p6.4 = 4.5, and the solid 
line is the best fit p 6A = 4.33 ±0.05, implying (3 = 0.94 ±0.02 
in the SL mode, very close to the value f3 = (7/9) 1/3 « 0.92 
suggested in ref. |}(| . For different pairs (p, q) the universality 
and the agreement with SL and LN is correspondingly good. 
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